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Abstract 

The density distributions of the two components of a trapped, ultracold Fermi gas with popu- 
lation imbalance reveal the effect of imbalance on superfluid pairing. We develop a path-integral 
derivation of the density, that takes into account both fluctuations beyond mean-field and effects of 
nonzero temperature. The calculated density profiles compare favorably to the measured density 
profiles, and illustrate the necessity to incorporate both quantum fluctuations and finite temper- 
ature effects. The temperature dependence of the density profiles, especially near the superfluid- 
normal interface, allows determining the temperature of the superfluid core in current experiments. 
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Recent progress in the trapping of cold atoms has attracted great interest related to a 
wide variety of fields: condensed matter physics, atomic, molecular and optical physics, 
astrophysics and physics of quark and nuclear matter. There exists a deep analogy between 
the dynamics of cold atoms, astrophysical systems, nuclear and plasma systems 0, 
3]- The possibility to tune the effective interaction strength and to control populations 
of different spin states provides a unique opportunity to investigate various phenomena of 
interacting many-body systems. Imbalanced fermi superfluids are in particular relevant 
for color superconductivity in dense quark matter [s], and for neutron-proton pairing in 



asymmetric nuclear matter 

In experiments jl, 0, [lO, 11, [H, a phase separation between the superfluid and 
normal component of an interacting Fermi gas of cold atoms with unequal spin populations 
has been observed. These experiments identify a shell structure of the fermion cloud, in 
which a superfluid core is surrounded with a Fermi gas in the normal state. 

The mean-field approach provides a convincing qualitative explanation of the phase sep- 
aration of imbalanced cold fermion atoms in a trap [3, 15, 16|. Mean field predicts a 



superfluid core (which at T = is completely unpolarized) surrounded with an imbalanced 
normal phase, and a discontinuous behavior of the majority and minority component den- 
sities at the phase boundary. The theoretical study of a trapped strongly interacting Fermi 
gas in the unitarity hmit at the zero temperature [13] shows a good agreement of the radii 
for the majority and minority components with the experimental results of Ref. [ll|. How- 
ever, by contrast with experiments, the boundary between the two phases provided by the 
mean-field approach (see Refs. [3, m, H 113 and the mean-field results below in the 



present work) is very sharp both for elongated and oblate traps, i. e., independently of a 
trap anisotropy, while experimentally observed phase boundaries are rather smoothed in the 
case of Ref. [10]. Here, we describe quantitatively this smoothing by taking into account 
effects beyond the mean-field approximation, at nonzero temperature. The density profiles 
derived in this work offer the possibility to determine the temperature of the superfluid by 
a fitting procedure, all other parameters being fixed by the experimental configuration. 

Recent experiments investigating pairing of cold atoms are realized in the crossover 
regime between weak coupling, where the Bardeen-Cooper-Schrieffer (BCS) theory holds, 
and strong coupling where the system is well described as a Bose- Einstein condensate (BEC) 
of molecules. In the BCS-BEC crossover regime, the critical temperature of the superfluid 
phase transition is substantially lowered with respect to obtained within the mean-field 
approach. The lowering of Tc is due to Gaussian fluctuations about the saddle point. The 
fluctuations influence also other parameters of the fermion system, such as the density. 
In addition to fluctuations, also other factors influence the density profiles calculated for 
cold Fermi gases: the possible violation of the local density a ppr oximation (LDA) in highly 



anisotropic traps observed by Hulet and co-workers [ill . |12| . [ISJ or the occurrence of non- 
trivial phases such as Fulde-Ferrell-Larkin-Ovchinnikov state as reported by Yoshida and 
Yip [18]. These mechanisms can be described using the method outlined below, but this is 
beyond the scope of the present treatment, which focuses on the role of fluctuations. High- 
resolution density distributions of cold ^Li atoms measured by Ketterle and co-workers [lOl] 
allow us to perform a quantitative comparison between the experiment and theory. 

There are different techniques to treat degenerate Fermi gases in the BCS-BEC crossover 
regime. The diagrammatic technique is applied by Perali et al. p!s|, Chen et al. [20|, and 
Taylor et al. |2l| to extend the results obtained by Nozieres and Schmitt-Rink (NSR) [23] 



for T = Tc to arbitrary temperatures. The NSR-like scheme (23| provides the equation of 
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state which is in an excellent agreement with Monte Carlo calculations within the whole 
range of the BCS-BEC crossover at finite temperatures. This approach is also capable to 
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show the universal thermodynamics of strong interacting fermions [2^ 

In order to investigate the phase separation of an interacting Fermi gas in a trap, 



yai 

use the path-integral formalism [26|, |27|, [28 



we 



The path-integral method is attractive because 
it allows to obtain reliable parameters for interacting fermions in the whole range of the 
coupling strength. However, it is necessary to go beyond the low-temperature approximation 



of Refs. [27|, |28| , because the contribution to the density due to fiuctuations is sensitive to 
the temperature variations, even in the low-temperature region. 

In the present work, the path-integral approach is extended to the case of arbitrary 
finite temperature and imbalanced spin populations, and applied to calculate fiuctuation 
contributions to the density. We consider a two-component Fermi gas with the spin states 
|-|-) and |— ) in an anisotropic parabolic trap within the LDA. The coordinate dependence 
of the density of fermions in a given spin state is determined through the chemical potential 



(1) 



characterized by the confinement frequencies uji and uj2. For given values of the local density 
n = -|- n„ and of the local density difference 5n = — n_, the averaged chemical po- 
tential yU = + fi_) /2, the imbalance potential ( = — /i-) /2 and the gap parameter 
A are found as a joint solution of the gap equation and of the two number equations. The 
thermodynamic potential and therefore the densities are determined in the quadratic ap- 
proximation with respect to the quantum fluctuations about the saddle point for arbitrary 
temperatures and taking into account the population imbalance. 

Within the NSR scheme, the Gaussian fluctuations do not feed back into the saddle-point 
equation for the functional integral, and the fluctuations only contribute to the fermion 
density and to the number equation. Hence the gap equation takes the same form as in the 
mean-fleld approximation. As shown in Ref. 2^, this is a natural approximation within 
the path-integral formalism. 

For a balanced gas, the saddle-point thermodynamic potential has a single minimum, 
which results in the gap equation. For an imbalanced gas, however, the saddle-point ther- 



modynamic potential can have two minima at A = and at A 7^ [29(]. Therefore, the 
phase boundary for an imbalanced Fermi gas in a trap cannot be determined from the gap 
equation, as well as from the Thouless criterion [s^], because they both allow one to flnd 
only a local minimum of the variational problem for the thermodynamic potential jlO]. As 
a natural extension of the NSR scheme to an imbalanced Fermi gas, we flnd coordinate- 
dependent values of the gap parameter through the straightforward minimization of the 
saddle-point thermodynamic potential [sH, 



sp 

v 



dk 

(2^ 



(3 



In (2 cosh PC + 2 cosh pE^) - 



|A|^ 
2P 
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(2) 



where a<j is the scattering length, ^i^ = — fi is the fermion energy, and = y + \A\ is 
the Bogolubov excitation energy. We use the units in which h = 1, m = 1/2, and the Fermi 
energy Ep = 1- 



3 



The fermion densities are a sum of mean-field and fiuctuation contributions. The current 
finite temperature implementation of the extended NSR scheme is similar to that used in 
studying d-wave pairing in Ref. j3l|. For the case of the s-wave pairing mechanism, we find 
that the fiuctuation contribution Ufi to the total density n is given by 



(iq 

(2^ 



im — ^, , . , au! 

g/3(<^+J7) _ 1 



(3) 



with (3 the inverse to the temperature. Here, no is an arbitrary positive integer, and the 
parameter 7 lies between two bosonic Matsubara frequencies flno < 7 < ^no+i, = 27m/ p. 
For the computation, a fast convergence of the integral over in ([3]) is achieved when 
7 = (2n + 1) tt//?, and an integer n is chosen in such a way that 7 ~ 1. The function J (q, z) 
of complex argument z is 



J (q, z) 



Mi,i (q, -z) '-^^ - Mi,2 (q, -z) '-^^^ 



Mi,i (q, z) Mi,i (q, -z) - (q, z) 



(4) 



The fiuctuation contribution 6nfi to the density difference 6n is given by the same expression 
as ([3]) with replacing the derivatives d/dfj, in the function J {q, z) by d/d(. The matrix 
elements Mj^^ (q? z) are given by 



Mi,i(q, z) 



X 



dk 

(2^ 
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2A;2 ' 2Ek (cosh /5Ek + cosh /5C) 
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Mi,2 {q,z) = -\A\^ J 



X 



(ik sinh /?i?k 
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1 

(2; + Ek - ^k+q) (^ + ^k + ^k+q) 



(6) 



where is the s-wave scattering length and -Ek = ^/^k + the energy of the Bogolubov 
excitation with 6]^ = k'^ — /x. Here, we use the units with the mass m = 1/2, ^ = 1, and the 
Fermi energy Ep = 1. These expressions incorporate not only the particle-pair and hole-pair 
excitations, but also particle-hole excitations, so that in the limit A ^ one obtains an 
interacting Fermi gas rather than the ideal gas. 
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The phase transitions of the imbalanced fermion 
for various thermodynamic variables (see Refs. j32l . 



;as were analyzed using phase diagrams 



The finite-temperature 



phase diagram pressure/temperature [37|, |38j is probably the first phase diagram which 



includes finite temperature and fluctuation effects as well as trap inhomogeneity. At a fixed 
scattering length a^, a temperature Tq of a tricritical point separates two types of the phase 
transition as follows. For T > Tq, the relative population imbalance 6n/n is continuous 
across the phase transition, so that the latter one is of the second order. On the contrary, at 
T < To, the relative population imbalance changes disco ntinuously at the phase transition, 
so that the latter one os of the first order [29}]. As a result, there exists the region where a 
uniform normal or superfiuid state cannot exist, and therefore in this region phase separation 



occurs 



The effect of fluctuations beyond mean-field on the phase diagram is not completely 
resolved [37|, |38|, although the effect of Gaussian fluctuations on the locus of the tricritical 



points was studied in [32|. In this contribution, rather than focusing on the hard task of 
correcting the mean-field phase diagram, we will investigate the effects of fluctuations on 
the density profiles, where the fluctuation contribution is straightforwardly obtained from 
Eq. iD. 

The phase separation exists at temperatures lower than the temperature of the tricritical 
point. Typical temperatures of the experiments on the phase separation of an imbalanced 
mixture of cold atoms are estimated to be as low as T ~ O.lTp, where Tp is the Fermi 
temperature Tp = Ep/kp. Therefore, the phase transition at the boundary of the superfiuid 
core in the experiments |10l . is expected to be of the first order. In order to get a definite 
answer to this question, we investigate the distributions of the gap parameter and of the 
density for a fermion gas in an anisotropic parabolic trap taking into account Gaussian 
fluctuations and using the parameters typical for the conditions in the experiment by Ketterle 
and co-workers ^IQ] . 

In Fig. [H (a), we plot ID profiles 6n obtained by the integration of 6n over y and z 
coordinates (corresponding to the profile of the integrated density difference from Fig. 3(d) 
of Ref. |10j). Fig. [!](&) shows the 3D distribution profiles along the lateral direction for the 



density difference, compared with the experimental distribution profiles from Ref. [10] . Both 
with and without fluctuations, there is a spatial phase separation, where a superfiuid core 
is surrounded by a normal phase. In the trap at a sufficiently low temperature the fermion 
system is almost exactly balanced in the core and almost completely polarized outside the 
core. Local densities for both majority and minority components change their values sharply 
when passing the phase boundary. Correspondingly, the total fermion density and the local 
density difference reveal discontinuities at the phase boundary. Those discontinuities are a 
consequence of the fact that the phase transition is of the first order. 

Within the BCS-BEC crossover regime, the Gaussian fluctuations lead to a substantial 
decrease of the core radius with respect to that obtained in the mean-field approximation. 
The fluctuations reduce the amplitude of the discontinuous change of the density difference 
and provide a partially smoothed (with a relatively small discontinuity) density profile near 
the phase boundary. A smoothing of the radial density profiles owing to fluctuations was 



also obtained by Chien et al. |37|, |38j. The profile for 5n, obtained taking into account fluc- 
tuations, lies drastically closer to the experimental radial profile [13] for the local population 
imbalance than the mean-field result. As well as the radial distributions, the ID profiles of 
the integrated density difference calculated taking into account Gaussian fluctuations are in 
a better agreement with the experiment ^] than the mean-field results in what regards the 



5 




FIG. 1: (color online), (a) The distribution profiles for the integrated density difference compared 
with the experimental data (black curve) from Fig. 3(d) of Ref. [13] . (b) The 3D distribution 
profiles along the lateral direction for the local density difference compared with the experimental 
3D distribution profile (black curve) from Fig. 3(e) of Ref. Dashed and dotted curves 

correspond to the present approach taking into account Gaussian fluctuations, for T = O.lTp and 
T = O.OGTp, respectively. Dot-dashed curves correspond to the mean-field approach. 



core size and the density profile. With respect to the present approach and to the exper- 
iment, the mean- field approximation overestimates the size of the superfluid core and the 
magnitude of the jump of the density difference at the phase boundary. 

The 3D distribution profiles along the lateral direction for the gap parameter calculated 
within the mean-field approximation and taking into account Gaussian fluctuations are plot- 
ted in Fig. [21 As seen from the comparison of the spatial distributions of the gap parameter 
obtained within the mean-field approach and within the NSR scheme, the fluctuations re- 
duce both the core radius and the values of the gap parameter inside the core. At the 
phase boundary between the superfluid and normal phases, the gap parameter discontinu- 
ously changes its magnitude. When taking fluctuations into account, the first-order phase 
transition between the superfluid and normal states occurs at a higher magnitude of the 
gap parameter than within the mean-field approach. This is a consequence of the fact that 
fluctuations lower the stability of the superfluid state. 

As found in Refs. [21I . [321 . [37I . [38| . the NSR scheme leads to difficulties near the 



unitarity limit, where the fluctuation contribution to the fermion density is not small and 
therefore can be hardly treated as a perturbation. For example, close to the critical temper- 
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FIG. 2: (color online) 3D distribution profiles along the lateral direction for the gap parameter 
calculated without (solid curve) and with (dashed curve) Gaussian fluctuations for cold fermions 
in the parabolic trap at T = O.lTp and 6N/N = 0.58. 



ature, the NSR approach predicts a bend-oyer behavior of the gap parameter as a function 
of the temperature 
matrix 
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showed that near unitarity, the susceptibihty 
\ 

(7) 

within the NSR scheme is not positive semi-definite. They interpreted this result as a 
breakdown of the NSR treatment. However, this breakdown occurs when the fluctuation 
contributions are sufficiently high - at temperatures close to Tc for a uniform Fermi gas. 
Also the bend-over behavior of the gap parameter predicted by Taylor et al. 2l|] occurs near 
the critical temperature for a balanced uniform Fermi gas Tc ~ 0.23Tp 26| . 

The fluctuation contribution to the fermion densit y st rongly falls down with decreasing 
temperature. The experiments of the MIT group (lol.]ll| are performed for lower tempera- 
tures T < O.lTp, at which the fluctuation contribution to the density is smaller than at Tc 
for a balanced uniform gas. In this connection, we can expect that the fluctuation contri- 
bution to the fermion density is relatively small with respect to the mean-field contribution, 
and therefore the extended NSR scheme can be applicable to the experimental conditions 
of Refs. jl, 0, [13, 11, 12, 13]. To verify this assumption, we study the behavior of the 
susceptibility matrix for cold fermions in a trap. 

The positive semi-definite matrix means that its eigenvalues are non-negative. The eigen- 
values Al and A2 of the matrix ([7]) can be expressed through the derivatives of the total 
density n and of the density difference Sn using as the independent variables the averaged 
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chemical potential /i and the imbalance potential (: 




djSn) 



+ 4 



dn 



(8) 



As follows from the relation 



djSn] 



(9) 



in the limit of a vanishing imbalance, when 6n = 0, 



0. Therefore in this case 



the eigenvalues of the susceptibility matrix are reduced to the derivatives ( |^ ) and 

fd{Sn)\ 

In order to check whether the susceptibility matrix is positive semi-definite at the con- 
ditions of the experiments under consideration we analyze the eigenvalues of the 
susceptibilit y m atrix as a function of the radius using the values of parameters relevant to the 
experiment [10|]. The eigenvalues are calculated taking into account Gaussian fluctuations 
about the saddle point. The 3D distribution profiles for the eigenvalues of the susceptibility 
matrix ([8]) are shown in Fig. [3l 




FIG. 3: (color online) 3D distribution profiles for the eigenvalues Ai, A2 of the susceptibility matrix 
for imbalanced interacting fermions in the parabolic trap at T = O.lTp and 5N/N = 0.58. 



The susceptibility matrix eigenvalues, as well as the gap parameter and the density pro- 
files, behave discontinuously at the phase boundary. We see that in the stable superfluid 
and normal states, both Ai and A2 are positive. This result shows that for the experimental 
conditions of Refs. [lo|, the NSR scheme extended to the imbalanced case provides a 
reasonable interpretation of the phase separation of a trapped Fermi gas. 
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Currently, the temperature in the experiment is determined by fitting to the tails of the 
remnant thermal distributions. One of the implicit assumptions is a good thermalization 
between the superfluid and the thermal cloud, and this assumption has recently been under 
further investigation j39|. Owing to fluctuations, the density profile is rather sensitive to 
temperature variations at low temperatures T <^ Tp. The good agreement between the 
temperature-dependent density profiles calculated here and the experimental density profiles 
allows us to estimate the temperature of the superfluid core. 

In conclusion, we have investigated the density profiles of a phase separated imbalanced 
Fermi gas in an anisotropic parabolic trap, taking into account Gaussian fluctuations near 
the saddle point at a nonzero temperature. Fluctuations lead to a smoothing of the density 
profiles and to a decrease of the radius of the superfluid core with respect to those predicted 
by the mean-field approximation. The calculated density profiles are in a good agreement 
with the experimental density profiles observed in Ref. [lO| and their temperature depen- 
dence enables temperature determination. The present approach allows us for a much more 
convincing interpretation of the experimental results on the phase separation of imbalanced 
fermions with respect to the mean-field approximation. 
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